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Abstract 

Defining a most economical parametrization of time-dependent B — ► p^ir^ decays, including a 
measurable phase a e g which equals the weak phase a in the limit of vanishing penguin amplitudes, 
we propose two ways for determining a in this processes. We explain the limitation of one method, 
assuming only that two relevant tree amplitudes factorize and that their relative strong phase, 
8t, is negligible. The other method, based on broken flavor SU(3), permits a determination of a 
in B° — > p^tt^ in an overconstrained system using also rate measurements of B 0,+ — > K *ir and 
B 0,+ — ► pK. Current data are shown to restrict two ratios of penguin and tree amplitudes, r±, to 
a narrow range around 0.2, and to imply an upper bound |a e g — a\ < 15°. Assuming that 8t is 
much smaller than 90°, we find a = (93 ± 16)° and (102 ± 20)° using BABAR and BELLE results 
for B(t) — > p^Tr^. Avoiding this assumption for completeness, we demonstrate the reduction 
of discrete ambiguities in a with increased statistics, and show that SU(3) breaking effects are 
effectively second order in r±. 
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I. INTRODUCTION 



A proposal made fourteen years ago to measure the CKM angle a in B — > -kit through an 
isospin analysis |l( was followed shortly afterwards by a suggestion to apply this technique 
also to the quasi two-body decays, B — > pir DD The analysis requires a construction of 
two pentagons, for B and B, the sides of which describe decay amplitudes into different- 
charge fm final states. This is a challenging task, requiring precise measurements of decay 
rates and asymmetries in all five modes. In addition it also involves a large number of 
discrete ambiguities and in certain cases continuous ambiguities in a j^j]. A simplification 
occurs when the decay amplitude of B -> A is 

much smaller than the amplitudes of 
the other four processes, in which case the pentagons turn into approximate quadrangles. 
Recently the BELLE collaboration reported evidence for B — > at a level implying 

that this simplification may not occur in practice. This seems to indicate that a useful 
measurement of a using the full isospin analysis may be impractical even with super-B- 
factory-like luminosities 0. 

A complementary and more promising way of learning a in these decays is based on 

n 

performing also a time- dependent Dalitz plot analysis of B° — > 7r + 7r _ 7r° |6[. One uses the 
interference between two p resonance bands to study the phase differences between distinct 
amplitudes contributing to the decay. This raises issues such as the precise shapes of the 
tails of the Breit- Wiener functions, and the effect of interference with other resonant and 
non-resonant contributions 0. A complete implementation of this method requires higher 
statistics than available today. 

An interesting and more modest question, which may already be studied now using 
time-dependent decay measurements of B°(B°) — > p ± ix T by the BABAR and BELLE P] 
collaborations, is what can be learned about the weak phase a using this limited information 
alone. Since these processes involve more hadronic parameters than measurable quantities, 
further assumptions are required to answer the question in a model-independent manner. 
Flavor SU(3) HQ, 

a symmetry less precise than isospin, provides a suitable framework 
for an answer. SU(3) symmetry relates B° — > p ± 7r T to processes of the type B — > K*7r 
and B — > pK 

H 

Allowing for certain SU(3) breaking effects, which may be justified on 
theoretical grounds and tested experimentally, improves the quality of such an analysis. A 
recent application of broken flavor SU(3) to the considerably simpler case of measuring a in 
B°(t) — ► 7r + 7r~ was studied in Hi- 

In the present paper we extend the SU(3) analysis of B° — > 7r + 7r~ - to study B°(t) — > p ±r K T . 
We also suggest an alternative approach to measure a in B — > p ± 7i T , which does not rely on 
flavor SU(3), but instead reduces the number of hadronic parameters by two when assuming 
that tree amplitudes factorize to a very good approximation. 
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Several earlier studies of a in B — > p ± n T have been performed. An application of flavor 
SU(3) to these processes has been carried out in 0, 13- We study a wider range of 
aspects, such as consequences of factorization of tree amplitudes, SU(3) breaking effects, 
bounds on ratios of penguin-to-tree amplitudes, and the nature of discrete ambiguities in 
values obtained for a. For completeness, we will present a proof of several bounds on two 
immeasurable quantities, af s — a 14, 15], one of which will turn out to be stronger than 



bounds obtained earlier. We will then define a measurable phase a e g which provides an 
approximate measure for a. 

An SU(3) relation between B° — > p ± n ± and B° — > K* ±r K T has also been discussed in Q|; 
however this work made no use of the interference between B°-B° mixing and B — > pn decay 
amplitudes which is a crucial input in our study. Implications of B°(t) — > p ± 7r =F for a global 
SU(3) fit to all charmless B decays into pairs of a vector and a pseudoscalar meson, B — > VP, 
were studied recently in Our study will focus on B — > p ± 7r =F and on their direct SU(3) 
counterparts. Our model- independent approach differs from other studies of B(t) — > p ±r K T 
involving a priori calculations of decay amplitudes and strong phases based on QCD and 
factorization 



18, 



The paper is organized as follows. Section II provides definitions of decay amplitudes and 
expressions for time-dependent decay rates in terms of a minimal set of parameters describing 
these amplitudes. Section III defines and discusses the use of a e g , a measurable that is equal 
to a in the absence of penguin amplitudes. Section IV considers a seemingly useful method of 
reducing the number of hadronic parameters by assuming approximate factorization of tree 
amplitudes, pointing out its intrinsic limitation. Section V draws SU(3) relations between 
B — ► p ±r K T and several processes of the type B — > K*n and B — ► pK. In Section VI 
we summarize the current experimental measurements of relevant rates and asymmetries, 
deriving numerical bounds on ratios of penguin and tree amplitudes in B° — + p ± n T and on 
the shift a e Q — a, obtaining a range of values for a. In Section VII we study the sensitivity to 
experimental errors of the flavor SU(3) method for determining a. This discussion involves 
certain discrete ambiguities, which will be discussed briefly in this Section, and will be dealt 
with in more detail in an Appendix. We conclude with a summary in Section VIII. 

II. AMPLITUDES AND TIME-DEPENDENT DECAY RATES IN B — > p ± 7r T 

We start by setting notations and conventions. B decay amplitudes, A + and A_, are 
denoted by the charge of the outgoing p, and corresponding 5° amplitudes into charge 
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conjugate states are denoted by A + and A_, respectively: 



A, 



A(B° 
A(B° 



p IT 



A_ = A{B° 



A(B 



A, 



P 71 

p + 7T 



(1) 



Each of the four amplitudes can be expressed in terms of two terms, a "tree" and a "penguin" 
amplitude, carrying specific CKM factors. We adopt the c-convention, in which the top- 
quark has been integrated out in the b — > d penguin transition and unitarity of the CKM 
matrix has been used to move a V* b V u d term into the tree amplitude. Absorbing absolute 
magnitudes of CKM factors in tree (t) and penguin (p) amplitudes, we write 



A± =e ir >t ± +p ± , 
A±=e- i n ± + p± . 



(2) 



While dependence on the weak phase 7 is displayed explicitly, strong phases are implicit in 
the definitions of complex amplitudes. We define three strong phase differences, 



5± = arg (p±/t±) , S t = arg (t_/t+) . 



(3) 



For convenience we also define ratios of penguin and tree amplitudes in the two processes 
and a ratio of the two tree amplitudes, 



p± 




/ 


t± 


n 







(4) 



Counting parameters, we find a total of eight, consisting of seven hadronic quantities 
\t±\i \p±\> ^ and the weak phase 7 or a. We will assume (3 to be given [l5 , 20 1 and use 7 = 
71 — (3 — a. In general, the amplitudes t + (p + ) and £_(p_) have different dynamical origins and 
are expected to involve different magnitudes and different strong phases. Amplitudes with 
subscripts +(— ) in Eq. (J2J describe transitions in which the final-state meson incorporating 
the spectator quark is a tt (p) (cf. Fig. This characterization was shown to be useful 
in the context of an SU(3) analysis of charmless B decays into a vector and a pseudoscalar 
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meson, B — > VP 



12l Il7|, w here and t-(p-) represent SU(3) amplitudes (denoted 

tp(pp) and ty(py) in jlj . v\). This broader framework will be used in our discussion below. 

Let us now consider measurables in time- dependent rates. Neglecting the width dif- 
ference in the B° system, and neglecting tiny effects of CP violation in B°-B~° mixing, 
time-dependent decay rates for initially -B decaying into p ± ir T are given by [2l| 

T(B°(t) -> p ± 7r T ) = e- vtl - (|A±| 2 + \A T \ 2 ) [1 + (C ± AC) cos Amt - (S ± AS) sin Ami] , 

(5) 

where 

C±AC= 1 -^ L^- , S±AS = i ■ 6 

\A±\ 2 + \A T \ 2 \A±\ 2 + \A T \ 2 v ; 

Here T and Am are the average B° width and mass difference, respectively. For initially 5° 
decays, the cos Amt and sin Amt in (jSJ) have opposite signs. 

Counting the number of independent measurables, we find a total of six, consisting of 
two CP violating quantities, S and C, two CP conserving measurables, AC and AS", and 
two rates, (r±) = |(|A-|-| 2 + |^4 T | 2 ). These two rates are related to the CP conserving charge 
averaged p ± n T combined decay rate, T p7T , and the overall CP violating asymmetry, A^p, 

r<- = (r + > + (r_) , a& = |±j ~ £=j , (7) 

implying 

(r±> = in- (1 ± ■ (8) 

Of particular interest are the two direct CP asymmetries between B°(B°) — > p + n~ and 
B°(B°) — > p _ 7T + decay rates, 



"^CP — l-T 19 , I . 19 ' -^CP — i-r 19 , I A 19 ' V > 



pn 

CP' 



1^4+ 


2 _ 


^4+ 


2 


K 


2 + 


^4+ 


2 





2 _ 




2 




2 + 
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These may be expressed in terms of three of the above measurables, C, AC and A; 



^ P (1 + AC) + C ^(l-AC)-C 
^ CP ~ l+A&C + AC ' CP ~ 1-^pC-AC ' 1 j 

The above observables, of which six are independent, can be expressed in terms of the 
eight parameters describing B — > p ± 7r =F in (0) and (j3J). This leads to rather lengthy expres- 
sions, which we do not fully display. Here we give the example of the overall CP asymmetry, 
Aq P , which does not depend on S t , 

jp* = 2 smjfi + a) (r+ sin 5 + - r 2 r_ sin g_) 

CP 1 + r\ + rf(l + — 2 cos(/3 + a)(r + cos 5 + + r 2 r_ cos<5_) 
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As demonstrated by this example which contains ratios of amplitudes rather than the am- 
plitudes themselves, it is useful to consider ratios of rates, thus trading two CP conserving 
measurables for two parameters \t + \ and |t_|. Our study is simplified by a judicial choice of 
the remaining four observables, such that they depend only on the six parameters \r±\, 5±, 
S t , and a without depending on r t . 

We now display a convenient (but not unique) choice for this minimal set of observables. 
Two of the observables are naturally the direct asymmetries, Aq F , which depend neither on 
on r t nor on St, 

± 2r± sin 5± sin(/3 + a) 



CP 1 + r| - 2r ± cos S± cos(/3 + a 
Instead of S and AS we define: 

1 



(12) 



S 



AS 



2-^/(1 + AC) 2 - C 2 



2-^/(1 + AC) 2 - C 2 



I Apn \ V 2 / 1 Apn \ V 2 

1 T S\ C P \ i I Q acW ^CP 



(S + AS) (- 



1 + AT x 1/2 



l CP 



AC* 



-(S-AS)i- 



1 _ AfW \ V 2 



,4 p7r 
-^cp 



(13) 



(14) 



Note that S and AS are CP violating and CP conserving, respectively, in complete analogy 
to S, AS. They are free of r t , and their dependence on other hadronic parameters is given 
by 

1 



S 



sin 2a — (r+ cos 5+ + r_ cos S-) sin(a — /3) — r + r_ sin 2/5 cos(5+ — 



cos<5f 



(r + sin <5 + — r_ sin <5_) sin(a — /3) + r + r_ sin 2(3 sin(5 + — SJ) 



sin S 



cos 2a — (r + cos S + + r_ cos 5_) cos(a — (3) + r + r_ cos 2/3 cos (5+ — S-) 



(15) 
sin Si 



+ 



(r+ sin <5 + — r_ sin <5_) cos(a — (3) — r + r_ cos 2/5 sin(<5 + — <5_) 



cos St 



(16) 



where 



^J— = J (1 - 2r+ cos(/3 + a + 5+) + r 2 ) (l - 2r_ cos(/3 + a + <L) + r 2 ) . (17) 



In our discussion in Section VII of determining a in a broken SU(3) analysis we will use 
this most economical parametrization of time- dependent measurements in B°(t) — > p^n 11 , 
given by the four measurables, A^ p , S and AS" in Eqs. ()12jh ([To]) and (fTtjj) in terms of the 
six parameters, r±, S±, S t and a. 
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III. THE USE OF a eG 



We follow the simpler case of B°(t) — > 7r + 7r~, where a contribution of a penguin amplitude 
modifies the value of a to a e fr = \ arg[e _2l/3 y4(i? — > 7r + 7r~)74*(.B — ► 7r + 7r~)], measured from 
the two coefficients of the sin Ami and cos Ami terms j^j. The isospin analysis [1] provides 



a way of determining the shift a e g — a. In B — > p^vr 111 we now define two corresponding 
quantities 0, 14, 1^ |. 

«£ = ^rg (e-^A ± A* ± ) . (18) 

Note that these phases do not occur in the time-dependent rates (0) and are unmeasurable 
in B — ► p ± 7r T alone. Instead, the observables 5 ± AS © involve two other related phases 
which can be measured in these decays, 

2a% ±6 = arg (e- 2ip A±At ) = arcsin [ - S ^ AS ^ ] ( 19 ) 
eS V TJ {^1 - (C T AC) 2 J 

where 

6 = arg (A*_A+) (20) 

is an unknown relative phase between the two decay amplitudes. Consequently, although a^ s 
and a~ s cannot be measured separately, their algebraic average is measurable. We therefore 
define: 



1 / + -\ 1 
"cfr = - (« eff + a eff J = - 



S + AS \ . I S-AS 

arcsin — , + arcsin 



y/i - (c + Acy J yV 1 - (c- ac) 2 

' "(21) 

The two shifts a eS — a are expected to increase with the magnitudes of the corresponding 
penguin amplitudes, \p±\. A relation between — a, \p±\, 7 and corresponding charge 
averaged rates and CP asymmetries is readily obtained using Eqs. (J2J) (a similar relation in 
the different t-convention was shown to hold in B — ► n + 7i~ 123111 . 



4|p ± | 2 sin 2 7 =(|A ± | 2 + |A ± | 2 ) 



1- (^ P ) 2 cos2(a e ± ff -a 



(22) 



The left-hand-side of (|22jl can be bounded using flavor SU(3) as shown in Section V. This 
implies lower bounds on cos2(a^j- — a) or upper bounds on \af s — a\ (see Eqs. (|36|) . (jSSj) 
below). Using (J2TJ) will then provide an upper bound on \a e s — a\. 

IV. ASSUMING FACTORIZATION OF TREE AMPLITUDES 

Since the number of parameters in B°(t) — > p ± 7i Zf exceeds the number of measurables 
by two, a certain input is required in order to determine a from these measurements. This 
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input is provided by an assumption that the two tree amplitudes t± factorize and that their 
relative strong phase vanishes in this approximation. Given that factorization was shown 



to hold to leading order in l/m& and a; s (m&) in a heavy quark QCD expansion |2J, |2f 
we will proceed under this assumption. Thus, neglecting for a moment a ratio of two form 
factors contributing to t- and t + [19], we take r t = to be given by the ratio of 

corresponding decay constants, 

r t ~ j- = 0.63 , (23) 
Jp 

where f w = 130.7 MeV, f p = 208 MeV. We note that a value r t = 0.68 was obtained 
in a global SU(3) fit to all B — > VP decays supporting both factorization of tree 
amplitudes and the assumption that B — > 7r and B — > p form factors do not differ much 
from one another. The absolute value of \t + \ obtained in the fit of Ref. 17[ agrees with 



\t + /t\ ~ f p /fn, where t, the tree amplitude in B° — > it + tt~, is obtained from a global SU(3) 
fit to B decays to two charmless pseudoscalars j^J. This also supports factorization of tree 
amplitudes and an assumption that the B to 7r form factor varies only slightly with q 2 . 

Factorization of tree amplitudes also implies that to a good approximation S t ~ 0. A very 
small phase, 5 t — (1 ± 3)°, supporting this assumption, was calculated in (Somewhat 



larger values around —20° were obtained in the global fit 17].) Taking r t to be given by (|23|) 
and assuming S t ~ reduces by two the number of parameters describing B — > p ±r K T , to 
become equal to the number of observables. Although this situation seems perfectly suitable 
for a direct determination of a, we wish to point out its limitation. 

As noted above, the four observables -A^p; & anc ^ depend on six parameters, r±, 5±, 5 t 
and a, one of which is assumed here to vanish approximately, St ~ 0. The overall CP 
asymmetry Aqp, given explicitly in ()lip. provides a fifth measurable, depending also on r t , 
which is assumed to be given by (|23p. While in principle this permits a determination of 
a, this can be seen to rely on terms quadratic in r±. Expanding Eqs. (fTTj) and ()12j) up to 
terms linear in r±, we find 

= -A+ P + r t M C P + 0{rl) . (24) 

That is, at this order the three observables are not independent when r t is given. As we will 
show in Section VI, one expects r± to be small, r± ~ 0.2, implying that a determination of 
a using these assumptions will be very difficult. One may turn things around, however, by 
using the linear relation to determine r t and thereby test factorization. 



V. CONSTRAINTS FROM FLAVOR SU(3) 

Another way of adding an input into the analysis of B — > p ± ir T is provided by assuming 
flavor SU(3), as we show now. In order to improve the precision of our analysis, we introduce 
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SU(3) breaking corrections in tree amplitudes. These amplitudes, which can be shown to 
factorize to leading order in l/m& and a s (rrib) 24, will be assumed to involve SU(3) 
breaking factors given by ratios of meson decay constants. Penguin amplitudes, for which 
factorization is not expected to hold j^ElL be assumed by default to obey exact SU(3). 
The effects of SU(3) breaking in penguin amplitudes will be discussed further in Section fvTTl 
Strangeness changing amplitudes describing B — > K*ir and B — > pK will be denoted by 
primed quantities. The SU(3) counterparts of t+, t- and p±, (J2J), are given by 0, 01 

,, _ fx* Vub V us , fx* t , 
+ ~~ ~T~ v* v + ~~ ~T At + > 

JP v ub v ud Jp 

.i _ Ik V* b V us fxT, 
- ~ Tv*v ~ T ' 
v*v 

P± = T^r r P± = -\ V± , (25) 

V cb V cd 

where 

A = - \ 2/ = 0.230, ^ = 1.04, ^ = 1.22. (26) 

1 — A / ^ Jp Jtt 

SU(3) amplitudes represented by exchange and annihilation contributions (contributing 
to AS = and AS = 1 decays, respectively) are l/m& suppressed relative to tree and 
penguin amplitudes [2^ and will be neglected. We also neglect very small color-suppressed 
electroweak penguin contributions. These approximations and the SU(3) breaking factors 
in (j25I) can be tested in B° -> K* + R-, B + -> K* + K° and in other B — > VP decays 0. 
Other tests, relating CP asymmetries in B — > p^ir^ and in strangeness changing decays, will 
be discussed in Section I VI II Under these assumptions one finds the following expressions 
for strangeness changing decay amplitudes 

A(B + -> K*\ + ) = - A~V 
A(B + -> p+^°) = - A" V 

and 



assumpt 

HQ 



(27) 



A(B° -> p-/^ + ) =^\t_e^ - A _1 p_ . 



(28) 



Denoting charge averaged decay rates by r(JB — > /) = [r(S —>•/)+ r(5 — > /)]/2, we 
now define the following ratios of charge averaged rates, 

^0 = A 2 £(i? ^- + vr-) ^ + ^ AT(i? + ^ *°7r+) 

+ r(fi°^p+7r-) ' + r(fi°^p+7r-) ' 1 1 

T(B° -f p-7T + ) " r(fi°^p-7T + ) 
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where superscripts and subscripts denote the charges of the B and p mesons. Using Eqs. flT}- 
(JHJ, ()27|) and (|28jl. the following expressions are obtained in terms of the hadronic parameters 
r± and 5± and the weak phase (3 + a: 

n o = rj + 2r ± \\z± + Aj ^ + = r| 
± 1 - 2r±z± + 4 ' ± 1 - 2r ± ^± + ' 1 J 

where 

z± = cos5±cos(/? + a) , A + = W^A = 0.235 , A_ = A = 0.254 . (32) 

yip y /t 

Each of these four measurables provides an additional constraint on appropriate pa- 
rameters. [CP asymmetries in B° — > i^* + 7r~ and 5° — ► p~K + do not provide additional 
information, but can be used to test SU(3); see Eqs. ()60|) and (joT|) below.] This leads to 
an overconstrained system from which a can be determined. That is, the six observables, 
given in Eqs. (JUJ), (|T5|l and (fTEJ) and in one pair of equations (pTTj). can be used to solve for 
the six unknowns, r±, 5±, 8 t and a, as discussed in more detail in Section VII. In the present 
section we study bounds on r± and on a e R — a which follow from the four observables TZ± + . 

Each of the four expressions (|3*TJl may be inverted to write r± in terms of z± and a 
corresponding ratio of rates, 



(TZ° ± + A±) 2 4 + (1 - ft° ± )(74 - A 4 ± ) - + A±) z± 

o 



7^ 



^± 2 4 + (i - ft±)ft± - ft± z± 



i-rc 



(33) 



± 



The four expressions are monotonically decreasing functions of z± having their minima and 
maxima at z± = 1 and z± = —1, respectively, 



, < r± < — == 

i + Jni i - 



05 <^<r^- (35) 



i + y/nl r l - 

Using current constraints on 7 [15], 38° < 7 < 80° (at 95% confidence level), the lowest 
and highest allowed value of z± are -0.79 and 0.79, respectively. This determines slightly 
smaller ranges of r± than given by (|34|) and ()35|) in terms of measured values of 7Z± + . 

We note that one may use ratios of separate rates for B or B mesons instead of the ratios 
of charge averaged rates defined in ()29|) and ([30)1 . The above considerations and the bounds 
on r± apply almost equally to these ratios. Instead of factors z± = — cos 5± cos 7 one now 
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has factors cos(<5± — 7) or cos(<5± + 7), which are constrained to lie in a range between —1 
to 1. These then imply bounds of the form (J53jl and (jHSJ). For given measurements of rates 
and asymmetries, as specified in the next section, one may then compare the three types of 
ranges obtained for r± and choose the most restrictive ones. 

The four strangeness changing processes (}2T)) and (J25|) . which are expected to be domi- 
nated by penguin amplitudes, can also be used to set an upper bound on \a e s — a\. For the 
two charged B decays one has 



cos2(a± 



a) 



l-27^sin 2 (/3±a) „ 1-27^ 



A±2 



A±2 



(36) 



For the two processes involving neutral B decays one finds 
X 2 T(B° -> K*+ir-] 



\P+\ 2 
X 2 T(B° -> p-K~ 

\P-\ 2 



1 + X%r~ 2 + 2X 2 + r~ 1 z + > sin 2 7 



1 + XtrZ + 2\ 2 _rZ z~ > sin 2 7 



(37) 
(38) 



where the two inequalities follow simply from the identity and the inequality 1 + x 2 + 
2x cos S cos 7 = 1 — cos 2 S cos 2 7 + (x + cos 5 cos 7) 2 > sin 2 7. Combining these inequalities 



with 



we find 



lill5] 



cos2(o4 



a) > 



2TZ° ± 



1 *A CP 



(39) 



Thus, measured branching ratios and asymmetries, appearing on the right-hand-side of (J36|) 
and and listed in the next section, provide upper bounds on \af s — a\ and, using (|2Tj). 
they imply upper bounds on |a e ff — a\. 

VI. CURRENT RATES, ASYMMETRIES AND BOUNDS ON r± AND a eS - a 



The current measured branching ratios and asymmetries in B° — > p^ir^ and in SU(3) 
related processes are summarized in Table I For ratios of B + and B° decay rates we will use 



the lifetime ratio |2jj t(B + )/t(B°) = 1.077 ±0.013. The BABAR 
collaborations measured in B°(t) — > p ± 7r T also the four quantities 
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291 and BELLE 



C 



0.34 ±0.12 
0.25 ±0.17 AC 
0.31 ±0.10 



0.15 ±0.11 
0.38 ±0.18 S 
0.21 ±0.10 



-0.10 ±0.15 
-0.28 ±0.25 AS 
-0.15 ±0.13 



0.22 ±0.15 

-0.30 ±0.26 

0.08 ±0.23 (S = 1.7), 
(40) 
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TABLE I: Experimental charge averaged branching ratios and CP asymmetries of selected AS = and AS = 1 
B meson decays. For each process, the first line gives the branching ratio in units of 10~ 6 , while the second line 
quotes the CP asymmetry. Note, that the averages for p + ir~ and p~ir + final states were obtained using also the 
value for the summed branching ratio from CLEO and are thus not a simple average of BaBar and Belle columns. 

Mode CLEO BABAR BELLE Avg. 

B° -> p ± ttT 27.6i?;! ± 4.2 [2Q] 22.6 ± 1.8 ± 2.2 [g] 29.1±|;^ ± 4.0 [^lj 24.0 ± 2.5 

-0.088 ± 0.049 ± 0.013 [29] -0.16 ± 0.10 ± 0.02 [9] -0.10 ±0.05 

p + ir- - 12.7 ±2.0 19.5 ±5.0 14.2 ±1.9 

-0.21 ±0.12 -0.02 ±0.16 -0.16 ±0.09 

p-Tc+ - 9.9 ±1.8 9.6 ±3.4 9.8 ±1.5 

-0.47 ±0.15 -0.53 ±0.30 -0.48 ±0.14 





16±f ± 2 [22] 


11.9 ±2.0 [2a, Ml 


i4.sj:ti±i:^;| 


12.7 ± 1.8 






-0.03 ±0.24 a [23., .34J 




0.06 ±0.20 


p-R + 


i6.o±5;| ± 2.8 [go] 


8.6 ±1.4 ±1.0 [34] 


i5.i±!:!±l:f±I:? [35] 


9.9 ±1.9 6 






0.131°;^ ±0.04 ±0.13 [34] 


n 99+O.22+O.O6 rocrl 
u - zz -0.23-0.02 Uiil 


0.17 ±0.15 


B+ _> K*°TT + 


7.6±H ± 1.6 [30] 


10.5 ± 2.0 ± 1.4 [2JJ 


9.83 ± 0.90±i:§2 [22] 


9.8 ± 1.2 


p+K° 


< 48 [29J 






< 48 



°The two measurements of the CP asymmetry entering this average have opposite signs, Acp = 0.23±0.18^o'o6 

Q and A C p = -0.25 ± 0.17 ± 0.02 ± 0.02 Q- The combined error includes a scaling factor S = 1.9. 
6 The error includes a scaling factor S — 1.2. 

where the first values were obtained by BABAR, the second by BELLE and the third are 
their averages. Statistical and systematic errors were added in quadrature, and a scaling 
factor S = 1.7 is used for the error on the averaged value of AS. Using these values and the 
definitions (fT^jl and (flljl. one finds 

_ , -0.11 ±0.13 _ f 0.21 ±0.14 

S={ AS={ , (41) 

-0.17±0.18 [-0.19 ±0.19 

for BABAR and BELLE respectively. In view of the difference between the values of AS 
measured by BABAR and BELLE, we will not only take their average in the discussion 
below but will also treat them separately. 

Let us now consider SU(3) bounds on the penguin pollution parameters r± and af s — a. 
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Using the definitions (ffilj) and (|3U|) and taking branching ratios from Table I, one obtains 
the following values: 



n° + = 0.048 ± o.oio , n% = 0.034 ±0.006 , n°_ = 0.053 ± 0.014 , (42) 

but only an upper bound on IZt. Applying (|3*3*|) and assuming Gaussian distributions, these 
values lead to allowed ranges for r±. The most stringent bounds on r + follow from 
which gives at 90% confidence level: 

0.14 (0.16) < r + < 0.25 (0.22) . (43) 

Values in parentheses are obtained by using the central value of The most stringent 
bounds on r_ are obtained by using the ratios T(B° — > p~ K + )/T(B° — > p~7r + ) and their 
charge conjugates, instead of relying on the ratio of the charge averaged rates. Using the 
branching ratios and asymmetries of Table I, one finds the following range at 90% confidence 
level, 

0.14 (0.21) < r_ < 0.34 (0.29) . (44) 

The bounds and (jBj) are expected to be modified by additional SU(3) breaking 
effects which were not included in the analysis. In (J25j) we assumed exact SU(3) for penguin 
amplitudes. Somewhat smaller values of r± are obtained if SU(3) breaking enhances p± 
relative to p±, as it would, for instance, by assuming factorization for these amplitudes. The 
above bounds are somewhat wider than and, as expected, consistent with values obtained 
in a global SU(3) fit to all B — > VP decays Q, r+ = 0.17 ± 0.02 and r_ = 0.29 ± 0.04, 
obtained when \p + \ and \p-\ were not assumed to be equal. (Somewhat smaller values, 
r_ = 0.25 ± 0.03, were obtained in the global fit when assuming p + = —p-, as proposed 
in 00 Values on the low side, r+ = O.lOl^ and r_ = 0.10±g;g§, were calculated in QCD 
factorization 

Assuming ^4^ p w 0, Eqs. and imply 

| sin(a+ -a) | < yfjQ , \ sin(a+ -a)\ < sin(/3+a)| , | sinfe -a)\ < ^RP_ . (45) 

Note that in (|39p A^ P is 3.3a away from zero; nonzero asymmetries would improve the above 
bounds. Currently we find the following upper limits at 90% confidence level, 

K% \at s -a\ < 14.2° , K% \a+ s -a\ < 7.3°-11.7° , TZ°_ => \a~ s -a\ < 15.4° , (46) 

where the central upper limit (obtained by using 38° < 7 < 80°) is shown to improve slightly 
as a increases in the range 78° < a < 122° [15]. Using (}2"T]) . the second and third upper 
limits imply 

|cn eff - a\ < 11.3° - 13.5° , (47) 
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where the bound is improved slightly as a becomes larger within the above range. A recent 
study of a in time-dependent asymmetries in B° — > 0] favors the upper part of this 

range. 

The upper bound (|47|L which may be improved by a few degrees through more precise 
measurements of branching ratios, including a first observation of B + — > p + K°, is expected 
to be modified by SU(3) breaking. For instance, if SU(3) breaking enhances p'± relative 
to p± (as it would by assuming factorization for these amplitudes), then the upper bound 
becomes stronger. In this respect these upper limits may be considered conservative. In any 
event, even if SU(3) breaking suppresses p'± relative to p± by 20 or 30 percent, one expects 
the upper bounds to change by this amount. This result provides an important conclusion, 
implying that in time- dependent decays B° — > p ± 7r =F a may be measured through a e g with a 
precision of about ±15°: 

|a c ff — a | < 15° (including SU(3) breaking) . (48) 

This accuracy is comparable to that of measuring a through a e g in time-dependent B° — > 
p + p~ decays L Here the shift caused by the penguin amplitude is constrained in the 
isospin analysis l|, |42j,|43j by measured branching ratios of B° — > p + p~ , B + — > p + p° and by 
an upper bound on B° — > p°p° [44J to a range, \a e g — a\ < 17°, at 90% confidence level. 

Using Eq. (JUj), the current data pU)l may now be translated into solutions for a e g. In 
order to reduce discrete ambiguities caused by the few branches of the two arcsin functions 
in (J2IJ), we note that 

(2a+ +5)- (2a c - ff -5) = 25 t + 0(r ± ) . (49) 

We will make a conservative assumption that the two angles on the left hand side differ by 
much less than 180°, which is equivalent to assuming that St is much smaller than 90°. (Note 
that QCD factorization predicts a very small value, St — (1 ± 3)° 3], while a global SU(3) 
fit finds S t — —20° 0| •) This mild assumption can be checked experimentally by measuring 
the phase difference a.rg(A_JA + ) using the overlap of the p + and p~ resonance bands in the 
B° — > 7T + 7r~7r° Dalitz plot 29]. (This measurement is expected to be feasible much before 
a complete isospin and Dalitz plot analysis [f| can be performed.) The measurable phase 
difference a.rg(A_ / A + ) is dominated by S t = arg(£_/£ + ). The difference, \axg(A_ / A + ) — S t \, 
is governed by the subdominant amplitudes p±. We have checked that this difference is 
less then 25° at 90% confidence level, when r + is in the range flUiJ), consistent with (p^ 
and (j4*2*|) . and when r_ is in the range (j44j) . consistent with (|3*3*j) where R°_ is replaced by 
T(B° — > p~ K + )/T(B° — > p~7i + ). Therefore, a small measured value of &rg(A_/A + ) would 
imply that St is much smaller than 90°, confirming our assumption. 
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Applying (|21|) separately to the BABAR and BELLE measurements, (which differ in their 
AS* values by 2a) and to their averages, we find that in all three cases the experimental errors 
in C, AC, S and AS translate into quite small errors in a e ff, ±5°, ±13° and ±4° respectively 
In each case one finds for the central values of these measurements only two solutions in the 
range < a e s < n: 

BABAR: 93°, 177° 
«eff = < BELLE : 102°, 168° • (50) 
Average : 94°, 175° 

Excluding by a + (3 < tc the three values near 180°, corresponding to an ambiguity a — ► 
37r/2 — a, we find 

(93 ± 5 ± 15)° BABAR 
a =< (102 ± 13 ± 15)° BELLE , (51) 
(94 ± 4 ± 15)° Average 

where the first error is experimental and the second is theoretical, coming from the bound 
(j4*Kj) . Note the weak dependence of a on AS, the central values of which have opposite signs 
in the BABAR and BELLE measurements (I40jl . Combining for simplicity the experimental 
and theoretical errors (|5T| in quadrature gives for the average 

a =(94 ±16)°. (52) 

All the above results are in goo d agreement with the range 78° < a < 122° obtained 
from other CKM constraints |15| . For comparison, we note that the world average CP 



asymmetries measured in B° — > tt + tt~ 45j, |46| have been recently studied in Ref. 13j and 
were shown to imply a comparable range, a = (103 ± 17)°, favoring large values of the weak 
phase in this range. 



VII. EXTRACTING a 



As mentioned, the observables in B —>■ p ± 7r =F and the SU(3) related rates are sufficient 
for determining a. One has four independent observables S, AS, A^ P in time-dependent 
decays B — > p ± vr T , and additional four independent observables 1Z° ± , 1Z± in AS = 1 decays 
(where currently only an upper bound on TVl exists). These eight observables provide 
an overdetermined set of conditions, as they depend on only six parameters, the hadronic 
parameters r ± , 5±, 5 t and the weak angle a. The set of Eqs. (H2), (fT3j). (fTBj) and (|3*T|) 
then allows for an extraction of a as well as all the hadronic parameters. No assumption is 
required about 5 t , thus relaxing the mild and experimentally testable assumption made in 
the previous section in order to obtain the unambiguous ranges (JS"Tj) . 
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FIG. 2: Confidence level (CL) as a function of a for a generated set of data using a choice of parameters, 
r+ = 0.18, r_ = 0.23, 6+ = 30j\ 5- = -55°, S t = 170° and a inpu t = 100°. The interpretation of a confidence 
level is the same as in Ref. when Gaussian errors are assumed. Errors used for x 2 are the currently 
measured ones [yellow (light gray) region], those anticipated with ten times statistics [cyan (gray)] and 
hundred times statistics [purple (dark gray)] . We assume an experimental error in as in 1Z°_ . 

A solution of Eqs. (fT2"j) . (fTKJ) . (fT7)|) and (jSIJ) under these general conditions is shown for 
illustration in Fig. El which plots confidence levels (CL) as functions of a for different levels 
of statistics. To obtain the plots we generated data for the observables 5*, AS, A^ P , 1Z° ± 
and 1Z±, using a particular choice of values for the parameters 5±, r±, 8% (as specified in the 
figure caption) and an input value a = 100°. The errors on the observables were taken to be 
the currently measured ones, apart from an error on 1Z_, for which the error was taken to be 
the same as the current error on 7Z°_. Improvements in confidence level are shown for ten and 
hundred times more data than available today. One sees that with enough statistics only 
one solution at a = 100° survives in the range 0° < a < 180°. That is, for this particular 
choice of parameters one hundred times the present statistics implies an uncertainty of ±2° 
in the single value of a extracted at 95 % CL. We checked that the situation presented in 
Fig. |21 is generic and applies to a large range of hadronic parameters. 

The ambiguities in a, which are eventually resolved with high enough statistics, are seen 
in Fig. 2 to imply a large range of allowed values of a at current statistics. To get a quick 
insight into the origin of these ambiguities, let us first explore the case of r± — 0. In this 
(oversimplified) case, involving merely mixing induced CP violation, the only observables 
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carrying information on a are S and AS, 



S = sin 2a cos St , 
AS = cos 2a sin 5f . 



(53) 



Assuming that S and AS are measured precisely, a solution for sin 2a is given by 
(sin2a) 2 = i jl + 5 2 - (AS) 2 ± ^/[l + £ 2 - (AS) 2 ] 2 - 4S 2 } . 



(54) 



The two signs in front of the square root of the discriminant correspond to the following 
map 

^ = {«-£-§> *t-f-2a}, (55) 
or equivalently to the following interchange in Eqs. (|53j) : 



sin 2a «-> cos <5 t , 
cos 2a «-> sin 5 4 . 



(56) 



The other discrete ambiguities of (J53)l are 



= {a -> a + tt} , (57) 



7T 



Pr/2 = -»• - - a, S t -> -<J t } , (58) 
P_ = { a _> - a , <J t tt - <J t } . (59) 

While none of these transformations relate to one another two values of a in the allowed 
range, 78° < a < 122° [15], combinations of these transformations, such as = {a ^ 

3n/2 — a, ...}, are relevant to this range. However, this ambiguity, as well as the others, is 
resolved once higher order terms in r± are taken into account, as can be seen in Fig. 2. (See 
Appendix A for details.) Namely, the complete set of equations, (|12|). (|T5j) . (|16|) and ()31j). 
is not invariant under these transformations, which are violated by terms of order r±. Note 
that although the ratios TZ± + are formally of order r± (because of the multiplicative factor 
A 2 in their definitions), they are in fact zeroth order. That is, they must be measured to an 
accuracy of order r± in order to resolve the ambiguities. It turns out that at least one pair 
of (R±, R±) ought to be measured to this precision. (See Appendix A.) 

Since the extraction of a relies on a given scheme of broken flavor SU(3), one may wonder 
how SU(3) breaking effects other than those included may affect the value of a. Flavor SU(3) 
is used to fix the values of r±, which we have shown to be small, r± ~ 0.2. The extracted 
values of a are given to zeroth order in r± by ([54)1 . Terms of order r± are affected by SU(3) 
breaking corrections, which by themselves are approximately r±. Therefore the overall SU(3) 
breaking effect in a is expected to be of order r\. 
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This is demonstrated in the following way. We first generate values for observables by 
randomly varying 5t, 5± and a, and taking r± to vary in the allowed ranges ()43j) and (|44j). 
In addition to SU(3) breaking in tree amplitudes we now allow also for the SU(3) breaking 
in penguin amplitudes by introducing positive parameters, 0.7 < c± < 1.3, and writing 
p'± = — c±\~ 1 p± instead of (|25[). [We neglect SU(3) breaking in 8± which is nonleading in 
the sense that the term 2r±\ 2 ±z± in the numerator of the first of Eq. ()31|) is smaller than the 
r\ term.] Generating data under this assumption, we then extract a using Eqs. (|12p. (|15|). 
(fTSJl and (jHU), where SU(3) was assumed to be present only in tree amplitudes. Running 
a Monte Carlo program for 10000 different configurations shows that the local minimum in 
X 2 shifts by only a/ ((a out — a [n ) 2 } ~ 2°, which is indeed of order r\. 

There exist experimental tests of flavor SU(3) and SU(3) breaking corrections, in terms 
of equalities between CP rate differences in SU(3) related processes. Two of these relations 
follow from Eqs. and (J2HJ) Q, 



T(B° 
T(B° 



p IT 



p TV 



r(5° 



t(b° 



p TV 



P TV 



A 

f K > 

u_ 



T(B° 



r(5° 



K* 



71 



-V(B° 
p + K~) - T(B° - 



p-R- 



(60) 
(61) 



(') ( f ) 

These relations test the equality of products \t± \ \p± \ sin 7 sin 5± in AS = and AS = 1 
decays. Using current values in TablelU one finds that, while the signs of the two asymmetries 
in the second equality confirm the SU(3) prediction, their absolute values differ by 2. Oct. 



VIII. CONCLUSION 

We have studied implications for a of time-dependent rate measurements in B — > p ± 7r T 
by proposing a parametrization which depends on a minimal number of hadronic parameters 
and observables. We have proposed one method based on factorization, which reduces by 
two the number of parameters to the number of observables. The limitation of this method 
was shown to be its sensitivity to terms quadratic in r±, the two ratios of penguin and tree 
amplitudes. 

Assuming broken flavor SU(3), which relates B — > p ± 7r =F to four processes of the form 
B°' + — ► K*7v and B 0,+ — ► pK, and using branching ratios measured for these processes, we 
calculated lower and upper bounds on r±, slightly below and slightly above 0.2. Defining a 
measurable quantity a e g, that becomes a in the limit of vanishing penguin amplitudes, we 
calculated upper bounds on |a e ff — <x\ in a range 11° — 13°, which are expected to be at most 
about 15° when including unaccounted SU(3) breaking effects. 

In order to resolve a discrete ambiguity in a, we assumed that the relative strong phase 
of two tree amplitudes, 5 t , is considerably smaller than 90°. This assumption, justified by 
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QCD factorization and by a global SU(3) fit to B — > VP decays, can be tested directly 
through a partial Dalitz plot analysis of B — > 7r + 7r~7r°. Using the BABAR and BELLE 
results for B(t) — > p ± 7r T this then implies single solutions, a = (93 ± 16)° and (102 ± 20)°, 
respectively, and an average a = (94 ± 16)°, taking into account an error scaling factor. 

Finally, using a complete set of measurables, including CP asymmetries and avoiding 
any assumption about 5t, we presented numerical studies demonstrating the feasibility of 
determining a and the reduction of discrete ambiguities with statistics. We have also shown 
that SU(3) breaking effects, which were not already included, are expected to be very small, 
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APPENDIX A: FURTHER DETAILS ON DISCRETE AMBIGUITIES 

Discrete ambiguities were found in Section VII when order r± corrections to observables 
were neglected. This lead to a 16-fold ambiguity on a in the range a G [0, 2tc) spanned 
by the transformations P t , P n , P n /2 and P_ given in ([57)1. (fo^j) and (|5^j). Let us now 
discuss how the higher order terms in r± affect the ambiguities. First of all, the (unphysical) 
transformation is an exact symmetry of Eqs. (fT2"|). (fT3j). (fTB^l . and (|3*Tj). if extended to a 
transformation on strong phases 



The other symmetry transformations Pi = {Pt, P-k/2, P-} receive higher order corrections. 
To see under which conditions they remain ambiguities, let us expand Eqs. (JUJ), (fT5j) . (JHJ), 
and (|3*T|) to first order in r±, where we count 8 t ~ r± ~ A <C 1, 

S — sin 2a = sin 2a [r + cos(/3 + a + 5+) + r_ cos(/9 + a + 

(A2) 

— sin(a — p) [r + cos 5 + + r_ cos 5_] , 



of order r\ 



P w = {a — > a + 7r, 5± — > 5± + 7r, 5 t 



(Al) 



A5 1 = cos 2a sin 5 t + (r+ sin 5 + — r_ sin 5_) cos(a — 0) 
Aq V = — 2r± sin 5± sin(/3 + a) , 



(A3) 
(A4) 




(A5) 
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The parameters r± are obtained from 



K% = sjTlt = r± . (A6) 

One is then left with six Eqs. (jA2|) . (|A3|) . (|A4|I and (|A5|) for four unknowns, 5±, 5 t , and a. 
In order to check for leftover ambiguities, let us assume that there exists at least one solution 
for (|A2|) . IjASj) . f)A4|) . and (|A5|I which we denote by 5°, and a . The transformations Pi 
may give us another viable solution for a, 

a aev > = P ia ° + 5a , (A7) 

together with new values for the other paremeters, 5^ w , £° ew . From (|A4|) we get (to leading 
order in the small parameters, r±, 5^ ew and 5a), 

sin 51™ = A Z P - - , (A8) 

± 2r± sin(/3 + P;a°) ' v ; 

while (|A5J) gives 

COSd± ~2r ± A|cos(/3 + P i a°) ' (AJj 
In general Eqs. (|A8jl . (JA9j) are not simultaneously satisfied which resolves the ambiguity. 

In case that 7Z± + is not measured to order r±, i.e. to a precision of about 20%, the 
ambiguity is retained if the right hand side of (jA8|) is not larger in magnitudes than one, 
leading to a solution for (jA2|) and (jA3|) : 



A ncw = 1 [AS - (r, sin 5" cw - r_ sin £* w ) cos^a - 0)1 , (A10) 

cos(2Pj« u i 



<5a 



— nx { S - sin 2a° + sin(P i a° - (3) \r+ cos <5" ew + r_ cos cT cw l 
2P,a u ) I 



2cos(2Piur*; v ( - A11 ^ ) 
- sin 2a° [r + cos(/3 + P^a + <^ ew ) + r_ cos({3 + P;a° + <T ew )] } . 

These expressions show the existence of further ambiguities in 5a of order r± caused by 
twofold solutions for 5^ ew in (jA8|) or (jA9|) . [Note that 5f™ and 5a are 0(r±) in accordance 
with our expansion.] This shows the importance of measuring TZ± + as precisely as possible. 
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